在分子动力学模拟中，有时需要固定一部分原子，例如，使用velocity命令对体系进行拉伸时，固定底部原子，对顶部原子施加一定的速度实现拉伸。

在纳米磨削中，也需要固定设置一个固定层，防止工件在磨削力的作用下移动。

在lammps中固定原子的方式有多种，我比较常用的是fix setforce命令。

fix setforce命令语法格式为：

其中：

fix setforce原理是在xyz方向上设置原子新的受力值，如果某个方向不需要设置，可设为NULL。

如果打算固定boundary组，就可以把boundary组各个方向的受力设为0：

如果只是使用setforce一句命令，有可能还会存在一个问题。

当boundary组的原子手里为0，但如果boundary组原子有非零温度，也就是原子有速度的话，即使受力为0，原子也会按照初始的速度继续移动，同样达不到固定原子的效果。

这种情况下，需要配合velocity使用，把该原子组的速度设为0.

上面两句代码配合使用，就可以固定boundary原子。

更多lammps教程和源代码请关注微信公众号：lammps加油站